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Within the framework of the current models for the slow neutron capture (s) process 
in Asymptotic Giant Branch (AGB) stars we explore the uncertainties introduced by the 
assumptions made on stellar modeling. On the basis of star models constructed with three 
different evolutionary codes we generate detailed neutron-capture nucleosynthesis post- 
processing models. The main difference among the codes is that one of them includes 
an overshooting mechanism. As a result, the neutron fluxes are stronger both during 
the interpulse periods and within thermal pulses. The efficiency of the 13 C source is also 
studied and we find that a linear relationship exists between the initial 12 C in the intershell 
and the maximum number of neutrons produced. 

1. INTRODUCTION 

The s process had been suggested in the 1950s to account for the nucleosynthesis of 
elements heavier than Fe The s-process path of neutron capture generally follows 
the valley of (3 stability, however, for some values of neutron density and temperature, 
branchings are open and nuclei outside the valley of stability can be produced. It can be 
shown that in equilibrium conditions the s-process distribution follows the rule (JaNa — 
const, where a a is the neutron capture cross section of the isotope A and Na its abundance 
in number. The rule is valid locally, far from neutron magic nuclei which act as bottlenecks 
of the flux and produce discontinuities in the distribution. An important parameter is 
the neutron exposure, r = Jq N D v t hdt, where v t h is the thermal velocity of the gas. The 
neutron exposure determines the final distribution of heavy elements; for higher neutron 
exposures heavier elements are produced 0. 

The major astrophysical site for the s process has been recognised to be the deep 
layers of AGB stars, which are low to intermediate mass stars in a late phase of their 
evolution. These stars experience recurrent convective instabilities in the region between 
the H- and He-burning shells (intershell) as a consequence of thermonuclear runaways 
of the He-burning shell (thermal pulses, TP). At the end of a TP, the convective enve- 
lope sinks into the intershell and dredges up material to the surface (third dredge up, 
TDU). A large amount of 22 Ne is present in the convective TPs as a product of the chain 
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14 N(a,7) 18 F( / 9 + z/) 18 (a,7) 22 Ne starting on the abundant 14 N from the H-burning ashes. 
However, the temperature in the TPs of low mass stars (M < 4 M ) is not high enough to 
significantly activate the 22 Ne(a,n) 25 Mg reaction. The 13 C neutron source is activated at 
lower temperatures (~ 0.8 x 10 8 K), however, an amount of 13 C higher than that present 
in the H-burning ashes is needed to reproduce the observed enhancements of heavy el- 
ements. To enable the formation of a region rich in 13 C, the 13 C pocket, some protons 
from the envelope must enter the 12 C-rich intershell. A favorable location for the occur- 
rence of the mixing is the sharp H/He discontinuity which is left over after TDU (see Fig. 
1), however, no standard stellar models have found the penetration of protons to occur. 
Recently, models including time- dependent overshoot |§, motivated by hydrodynamical 
simulations (e.g. |3J]), and models with rotation [H] have been able to create a proton-rich 
layer at the top of the intershell, after the end of the TDU. After less than a few thousand 
years 13 C and 14 N are formed in the region. Before the end of the interpulse period all 
13 C burns (a,n) and the region becomes enriched in s-processed material. At the end of 
the interpulse the pocket is engulfed by the following TP and is thus mixed with mate- 
rial from the previous TP and the ashes from the H-burning shell. Inside the TP, if the 
temperature at the bottom of the shell reaches 2.5 xlO 8 K, the 22 Ne(a,n) 25 Mg reaction is 
marginally activated and a second neutron flux can occur. After the TP the s-process-rich 
material from the intershell is dredged up to the envelope by the following TDU. This 
cycle is repeated over all TPs with TDU and the composition of the envelope throughout 
the AGB phase is changed by TDU and mass loss effects. 

Predictions from s-process models have to be compared with observational constraints 
from spectroscopical observations of chemically peculiar stars (see e.g. and with 

laboratory measurements of isotopic anomalies in silicon carbide (SiC) grains recovered 
from meteorites, and believed to be originated in AGB stars ||. Stellar models, together 
with Galactic Chemical Evolution models, should also reproduce the solar distribution of 
s-process elements [pi. For an extensive review on the s process in AGB stars see [110 . 



2. METHOD AND MODELS 



The s process has been computed making use of the Torino s-process code (TOSP, ||TT| ). 
This post-processing code calculates neutron captures on all nuclei up to Pb through 
the AGB phase. Since 13 C burns in radiative conditions, the 13 C pocket is simulated 
by specifying independent mass zones of varied extent and their 13 C amounts. The 
same 13 C pocket is adopted for all interpulse periods. Reaction rates for neutron pro- 



lyl for the C(a,n) O reaction, and from for the 



duction have been taken from 
22 Ne(a,n) 25 Mg reaction, excluding the contribution of the elusive resonance at 633 keV. 
In the temperature range of interest, from ~ 2 — 4 x 10 8 K, this rate is within a few 



percent of the one adopted in the NACRE compilation |TJ]]. Neutron capture rates are 



updated to the latest estimates (15 



The post-processing models have been generated for a star of 3 M and solar metallicity 
with inputs such as temperatures, densities, extent of convective zones and composition 
of the intershell taken from three different stellar evolutionary codes: the Italian Frascati 
RAphson Newton Evolutionary Code (FRANEC, |16|]), the Australian Mount Stromlo 
Stellar Structure Program (MSSSP, jl7fl ) and the German EVOL code Hl8| . Since the 
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Figure 1. Structure evolution thought two TPs of the 3 M© star model of solar metallicity as computed 
with the MSSSP and represented by the post-processing MOSN code (see §2). Convective regions are 
white and it is possible to recognise the convective zones of the 19 th and 20 th TPs, as well as the bottom of 
the convective envelope. The lines in the radiative regions represent the mass shells, which evolve in time. 
The region with the higher mesh density is the H-burning shell. The arrows give a good representation of 
the fact that the 13 C pocket is a local phenomenon involving a mass region of the order of half a tickmark 
on the y-axis, or less. The bottom of the convective envelope reaches 0.6643 M© during the 10 th TDU 
and the test exercises described in §2 have been performed in the region between 0.6638 and 0.6643 M©. 



codes are independent from one another, they contain several differences, e.g. with respect 
to the numerical treatment, opacities and nuclear reaction rates. The most important 
difference to be noted here, however, is the treatment of convective instabilities and 
the associated mixing. While the FRANEC code uses the Schwarzschild criterion to 
determine convective boundaries, the MSSSP code uses a special numerical scheme which 
may involve mixing of an additional stable mass shell during structure iterations. The 
EVOL code employs a time-dependent overshooting mechanism which leads to a deep 
penetration of the intershell convection zone into the C/O core during the TP (intershell 
dredge up, IDU), and to the formation of a 13 C pocket over a range of — 10 -5 M & . 
These effects depend on the value of the overshoot parameter /. The amount of TDU 
increases when moving from the FRANEC to the MSSSP to the EVOL code. 

In order to study the neutron flux in the 13 C pocket we have made use of a code 
specifically designed for light-element nucleosynthesis in AGB stars: the Monash Stellar 
Nucleosynthesis (MOSN) code, which uses the outputs of the MSSSP code. The nuclear 
network contains 74 light element species, up to the iron group, and 506 reactions updated 
as in [19| . Neutron captures on the missing elements are modeled by a neutron sink and a 
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fictional particle, g, is added to count the number of neutron captures occurring beyond 
62 Ni by an invented decay, with A=l s _1 : 62 Ni — > 61 Ni + g. Since the TOSP code does 
not includes proton capture reactions, the MOSN code has been useful in giving a good 
description of the neutron flux, particularly in regions of the pocket where the neutron 
poison reaction 14 N(n,p) 14 C and the consequent proton recycle effect is of importance. 
We experimented with the MOSN code, artificially introducing some protons below the 
H/He discontinuity left by the TDU to form a 13 C pocket. Then we let all the 13 C burn 
and checked the total neutron exposure. Since the maximum mass-shell resolution in the 
region of the pocket at the end of the interpulse period with the MOSN code is not smaller 
than ~5x 10~ 4 M , as shown in Fig. 1, to perform each test we put a constant amount 
of protons, rather than a profile, over a mass of 5 x 10~ 4 M Q . 

For the s-process calculations based on the FRANEC and MSSSP stellar models we 
used the same 13 C pocket profile as [|TI[] (see their Fig. 1), which has a total mass of 
5 x 10~ 4 M & and had been chosen in order to match observational constrains. The 
13 C pocket found self-consistently with the EVOL code and used in our calculation has 
been computed for the 5 th interpulse period with /=0.128 and extents over 1.7 xl0~ 5 
M Q . To compute the resulting neutron flux we performed the same tests with the MOSN 
code, however, this time we artificially introduced the 12 C, 13 C and 14 N found in the 
pocket with the EVOL code. We therefore specified the 13 C profile in the TOSP code to 
reproduce the r calculated with MOSN. 



3. RESULTS AND DISCUSSION 

3.1. The two neutron bursts 

The resulting 12 C, 13 C and 14 N profiles in the pocket as a function of the initial X p 
are presented in Fig. 2, as computed with the MOSN code. The total neutron exposure 
profile obtained when all 13 C has burnt is also plotted. Since the 14 N(n,p) 14 C reaction 
has a relatively high cross section ||20|| , nuclei of 14 N are a strong neutron poison during 
the s process. Some of the protons produced make more 13 C when recaptured by the 
12 C (recycle effect), others destroy the 13 C producing more 14 N. As a consequence, the 
total neutron exposure r grows, together with the 13 C profile, up to a maximum of 0.43 
mbarn -1 , and then decreases when the initial abundance of 14 N nuclei is comparable or 
higher than that of 13 C nuclei. This result confirms that the maximum value of r in the 
pocket chosen by []TTJ (see their Fig. 3 and 6) is typical for stars of solar metallicity when 
X( 12 C) ~ 0.2 in the intershell. We also made some tests changing the initial amount of 
12 C, and for different interpulse periods. As shown in Fig. 3 a nearly linear relationship 
exists between the final maximum value of r in the pocket and the initial 12 C, because 
the latter determines the 13 C and 14 N profiles and how many protons will be recycled. 
We note that the 12 C in the intershell varies from pulse to pulse. This linear relationship 
suffers from some uncertainties. While the temperature at the H/He interface at the end 
of TDU is almost constant (~ 10 7 K) for different interpulse periods, the temperature 
gradient towards deeper layers is much steeper for advanced TDU. The temperature at 
which 13 C and 14 N are formed in our calculations, at a mass ~ 3 x 10" 4 M below the 
discontinuity, is ~ 1.5, 2.2 and 2.5 xlO 7 K after the 1 st , the 10 th and the 17 th TDU 
respectively. The 12 C, 13 C and 14 N profiles can be slightly different and, as shown in Fig. 
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Figure 2. Resulting 12 C, 13 C and 14 N profiles in the 
pocket as a function of the initial number of protons 
introduced below the H/He discontinuity, after the 
end of the 10 th TDU (see Fig. 1). Also plotted is 
the total neutron exposure profile as computed at 
the end of the interpulse period after all the 13 C has 
burnt. 



Fi gure 3. Maximum neutron exposure in the 
pocket as a function of the initial 12 C. Squares, 
hexagons and triangles are results obtained for 
the interpulse periods following the 1 st , 10 th and 
17 th TDU respectively. Full symbols indicate the 
X( 12 C) value as computed by the MSSSP at such 
interpulscs. 



3, the T max versus 12 C relationship shows some spread. The smaller the extent of the 
pocket, the less the change of temperature gradient will effect the results. Future work 
includes more tests with the MOSN code, and s-process calculations with the TOSP code, 
using a different r max for different interpulse periods. 

As for the EVOL code, because of IDU, the initial mass fraction of 12 C is 0.48, and the 
resulting abundance profiles produce a r max of 0.88 mbarn -1 . The discrepancy with the 
maximum value of r we found for the same initial X( 12 C) with the test profiles, computed 
with the MOSN code, is possibly due to a slightly lower temperature at which the pocket 
has been computed, because of the smaller extent in mass, and mostly to the different 
choice of the 13 C(p,7) 14 N reaction rate. In the MOSN computation the latest rate by |2Tj 
has been adopted, rather than the previous [P2]j , which is 30% slower. 

To compare the neutron fluxes in the interpulse period and in the convective pulse, we 
plot in Fig. 4 the temperatures and neutron densities for the two neutron bursts, as com- 
puted by the TOSP code on the basis of the three different stellar evolution calculations. 
This comparison is only of qualitative character because we did not match the same TP 
number, core mass and envelope mass simultaneously for the three codes. However, the 
main finding is that the features of 13 C burning are fairly comparable in all cases (left 
panels in Fig. 4). All 13 C burns before the onset of the next TP because the temperature 
reaches 10 8 K (panel a). This is a very typical situation, however, it has to be noted and 
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Figure 4. The two neutron bursts as described with the different codes: temperature and neutron 
density as computed by the TOSP code in one of the zones that forms the 13 C pocket for a typical 
interpulse period (panel a and c) and at the bottom of a typical TP (panels b and d). The zero points 
in time represent: in panels a and c the time from the start of the interpulse period (about 10,000 yr) 
when T = 0.8 x 10 8 K; in panels b and d the time when T-^ ottom = 2.5 x 10 s K. Note the large differences 
in time scale, total r and neutron-density peak between the two neutron bursts. 



taken into account in future work that this might not occur during the earliest interpulse 



periods (see e.g. Fig. 4 of [18|). While the total neutron exposure is determined only 
by the initial amount of 13 C in any radiative layer of the pocket, the maximum neutron 
density can vary (panel c). However, it is so low that no branchings are activated. In the 
temperature range of interest the 13 C(a,n) 16 rate we use is about 50—70% lower than 



the one adopted in [FJ], however, this should not make much difference in the overall 
result. If 13 C burns faster the maximum neutron density will be higher, yet still unim- 
portant. In the TPs, the temperature profile has a narrow peak (panel b), and so has 
the neutron density (panel d). The neutron exposure is unimportant compared to that 
in the radiative layers, however, branching points are extremely sensitive to this neutron 
flux. Main differences show up between the FRANEC/MSSSP and the EVOL case. In 
the EVOL code, because of IDU, the temperature, and consequently the neutron density, 
stay high for a longer time and branchings are more active. 

3.2. s-process results 

We computed s-process models with the TOSP code over a total of 25 TPs with TDU 
for FRANEC, 18 for MSSSP and 13 for EVOL. In Table 1 the results for the envelope 
enrichment of s-process elements are compared with spectroscopic observations of AGB 
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stars and their relatives. We consider the parameters 'Is', which is the average of light 
s-process elements (Y and Zr) and 'hs', which is the average of heavy s-process elements 
(Ba, La, Nd and Sm). The ratios of these parameters with Fe, together with the ratios 
[hs/ls] and [Pb/Ba] are good indicators of the overall distribution of heavy elements, which 
is determined by the main neutron exposure from the 13 C pocket. Predictions are given 
for the last TP, initial values are solar (i.e. zero). To compare with measurements in SiC 
grains we consider three isotopic ratios: 138 Ba/ 136 Ba, 88 Sr/ 86 Sr and 96 Zr/ 94 Zr. The first 
two involve 138 Ba and 88 Sr, which are neutron magic nuclei, and are thus sensitive to the 
neutron exposure in the 13 C pocket, while the third ratio involves 96 Zr, which is produced 
through a branching at 95 Zr, thus it is sensitive to the neutron density in the TPs. The 
ranges of prediction are given for the last phase of the evolution when the condition for 
SiC grains to form, C/0>1, is satisfied in the envelope (last 5 TPs for the FRANEC, 4 
for MSSSP and 7 for EVOL). 



Table 1 

Results from the three computations as compared with some observational constraints, 

parameter measured FRANEC MSSSP EVOL 

[ls/Fe] a 0. -> 1.5 6 0M 095 (L52 

[hs/Fe] a 0. -» 1.5 fe 0.55 0.52 0.84 

[hs/ls] a -0.8 -> 0.5 b -0.40 -0.43 0.35 

[Pb/Ba] a - -0.52 -0.53 0.06 

5( 138 Ba/ 136 Ba) c ~ -250 d -289 -> -272 -331 -> -332 398 -> 647 

<5( 88 Sr/ 86 Sr) c 51±103 e 208 -> 294 20 -> 149 158 -» -293 

<5( 96 Zr/ 94 Zr) c -686 ± 187 / -505 -> -372 -718 -> -610 46 -> 1120 
a [X/Y]=log[(X/Y) meas /(X/Y) Q ]. b Range of values measured in stars of Z ~ Z Q from 0, typical 2a 

error ±0.25. c (5(x/y) = [(x/y) meas /(x/y) - 1] x 1000. d Bulk SiC grains from ||. e Average of single 
grains from [23]. Average of single grains from p5|. 



Since the results from the FRANEC and MSSSP computations, which use the same 
13 C pocket, are very similar, we can conclude that the important uncertainty in s-process 
predictions is the 13 C pocket, not the intrinsic differences among different codes (with 
no overshooting mechanism included). In the EVOL code, the neutron exposure in the 
pocket is higher and the production of the heavier s-process elements, Pb and the magic 
138 Ba are favored. The s-process element enhancements with respect to Fe, of about an 
order of magnitude in AGB star envelopes, are reproduced by the models. However, due 
to the uncertainties of the spectroscopic measurements on the observational side, and of 
the 13 C pocket on the theoretical side, we cannot draw strong conclusions about which of 
the models match the observations better. The ratio 138 Ba/ 136 Ba measured in bulk SiC 
grains indicates that, on average, the maximum neutron exposure in AGB stars should 
not exceed ~ 0.4 rnbarn . The ratio 96 Zr/ 94 Zr is matched when the neutron density in 
the TPs is not higher than about 5 x 10 8 n/cm 3 for a long time, a condition not satisfied 
by the EVOL model (see Fig. 4), which produces 96 Zr in a significant amount. 
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The same MSSSP s-process calculation was performed with a different 13 C profile, giving 
more of the total mass of the pocket to the zone with higher r (~ 0.4 mbarn -1 ) . Previously, 
as in ||11| , greater mass was given to the region with lower r (~ 0.15 mbarn -1 ). We found 
that [ls/Fe] and [hs/Fe] reached 1.27 and 1.0 respectively, while [hs/ls] and [Pb/Ba] kept 
to -0.27 and -0.48. The <5( 138 Ba/ 136 Ba) did not change much: -382 -> -368. The 
<5( 88 Sr/ 86 Sr), however, increased to 447 — > 651, going outside the measured range and 
supporting the profile weighted on the lower r. 
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